Quantum critical scaling and the Gross-Neveu model 
in 2 + 1 dimensions 
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Abstract - The quantum critical behavior of the 2 + 1 dimensional Gross-Neveu model in the 
vicinity of its zero temperature critical point is considered. The model is known to be renor- 
malisable in the large A*' limit, which offers the possibility to obtain expressions for various ther- 
modynamic functions in closed form. We have used the concept of finite-size scaling to extract 
information about the leading temperature behavior of the free energy and the mass term, defined 
by the fermionic condensate and determined the crossover lines in the coupling (g) - temperature 
(T) plane. These are given by T ~ |g — gc|, where gc denotes the critical coupling at zero tempera- 
ture. According to our analysis no spontaneous symmetry breaking survives at finite temperature. 
We have found that the leading temperature behavior of the fermionic condensate is proportional 
to the temperature with the critical amplitude -^tt. The scaling function of the singular part of 
the free energy is found to exhibit a maximum at ^ corresponding to one of the crossover lines. 
The critical amplitude of the singular part of the free energy is given by the universal number 
I [271^(3) ~ (f)] ~ —0.274543..., where (^{z) and Ch{z) are the Riemann zeta and Clausen's 
functions, respectively. Interpreted in terms the thermodynamic Casimir effect, this result implies 
an attractive Casimir "force" . This study is expected to be useful in shedding light on a broader 
class of four fermionic models. 



Introduction. — Quantum phase transitions take 
place at zero temperature by varying some nonthermal 
parameter, say g, such as composition or pressure, and 
are driven by genuine quantum fluctuations. At rather 
small (when compared to characteristic excitations in the 
system) temperature the singularities of the thermody- 
namic quantities are altered. It is then expected that the 
leading T dependence of all thermodynamic observables is 
specified by the properties of the zero-temperature critical 
points, which take place in quantum systems. In the close 
vicinity of a second order quantum phase transition the 
coupling of statics and dynamics introduces an effective 
dimensionality which depends upon (imaginary) time in 
addition to space [1] . In this case the inverse temperature 
acts as a finite size in the imaginary time direction for the 
quantum system at its critical point. This allows the in- 
vestigation of scaling laws for quantum systems near the 
quantum critical point in terms of the theory of finite-size 
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scaling [2-4]. 

The so called critical or thermodynamic Casimir effect, 
announced by Fisher and de Gennes in their study on crit- 
ical binary liquids [5] , is tightly related to the fluctuations 
of the order parameter in confined systems exhibiting tem- 
perature driven second order phase transitions in the bulk 
[4,6,7]. The named fluctuations give rise to a stress inside 
the confining walls, which impose effective boundary con- 
ditions on the system, depending on the behavior of the 
order parameter at the walls. The resulting stress can be 
described by the so called Casimir force, obtained as the 
derivative of the free energy with respect to the separation 
between the bounding walls. The critical Casimir force is 
a universal quantity, in the sense that it is independent 
of the details of the system. These are screened by crit- 
ical fluctuations that are correlated over a long distance, 
namely the correlation length, that grows indefinitely as 
we approach the critical point. By now it is well estab- 
lished that the Casimir force is attractive for symmet- 
ric boundary conditions (periodic, Dirichlet, Neumann), 
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while it is repulsive for non symmetric ones (antiperiodic, 
Dirichlet-Neumann) [4,6,8]. 

In analogy to the critical Casimir effect due to spatial 
confinement, we consider the critical behavior of quantum 
systems in the context of the critical Casimir effect, defin- 
ing the temporal Casimir effect [9,10] also dubbed Casimir 
effect in time [II], where the confinement is caused by the 
boundary conditions along the finite imaginary time di- 
rection. The singular part of the free energy of a d dimen- 
sional system in the vicinity of its quantum critical point, 
with hyperscaling satisfied, is expressed as 



(1) 



where g and v measure, respectively, the distance from 
the quantum critical point gc and the divergence of the 
correlation length. The number c is a non universal quan- 
tity encoding the irrelevant details of the quantum sys- 
tem. The function ^{x) is a universal scaling function, 
whose expression depends upon the universality class of 
the quantum critical behavior. At g — gc, the universal 
critical amplitude 2}(0) is the Casimir amplitude. Its sign 
indicates whether the fluctuation induced "force" is at- 
tractive or repulsive. In expression (I) we assumed that 
the dynamic critical exponent z — I. Let us mention 
that the critical Casimir amplitude in quantum systems 
is related to the specific heat and can be obtained from 
experimental data [4,12]. 

In this letter we tackle the quantum critical scaling be- 
havior in interacting massless fermions in order to make 
contact with the thermodynamic Casimir effect in time in 
such systems. To achieve our goal we employ the 0{2N) 
Gross-Neveu model [13] in 2 + 1 dimensions. The model 
describes discrete chiral symmetry breakdown. Its Hamil- 
tonian reads 



'lip 7 



d 
dxi 



■7 



_d_ 

dx-2 



where the fields i/> = xp{x,t) and 



less Dirac fermions with xpxp 



(2) 

= xjj{x,t) are mass- 

V'qV'qj 1^ are Dirac 
matrices and g - a. nonzero coupling between fermions. 

Model (2) may have applications that can be encoun- 
tered in a variety of physical phenomena spanning con- 
densed matter (superconductivity, cold matter, antiferro- 
magnetism) and elementary particles. The model is known 
to be exactly solvable in the large N limit, using varia- 
tional calculations [14,15]. It can be useful in shedding 
light on phenomena one might expect in more realistic 
models. 

The thermodynamics of (2) in the large TV limit and at 
finite temperature has been investigated in many occasions 
(See references [14, 15] and references therein). At zero 
temperature it exhibits a spontaneous symmetry break- 
ing at a critical coupling gc- Further, investigations at 
finite temperature in the presence of quantum effects have 
shown that, under certain assumptions related to the vari- 
ational problem (i.e., if only nonnegative solutions of the 



corresponding self-consistency or gap equation have been 
taken into account), a thermal driven phase transitions 
takes place at Tc > 0. In the neighborhood of Tc, the crit- 
ical behavior of a number of thermodynamics quantities 
has been determined. 

The aim of the present study is to go beyond these works 
and explore the mutual influence of quantum and classical 
fluctuations in the 2-1-1 Gross-Neveu model in the large 
N limit in the vicinity of gc keeping the temperature fi- 
nite (T — )• 0"'"). To this end we determine the scaling 
behavior of the correlation length and the free energy. In 
our analysis we release the constraint related to the vari- 
ational problem mentioned above. We will see that this 
leads to a "bona fide" solution for the free energy with 
reasonable physical output. To the best of our knowledge 
this is the first time ever that such a study of the model 
under consideration is performed. 

The remainder of this letter is organised as follows: In 
the next Section, we recall the main characteristics of the 
Gross-Neveu model and derive its critical properties in the 
vicinity of the quantum critical point (zero temperature 
critical coupling), keeping the temperature finite. To this 
end we compute the free energy and the gap equation to 
the leading order in the large N limit and analyse their 
corresponding behaviors as a function of the temperature. 
In the following Section we investigate the critical Casimir 
effect in time, by deriving the scaling function of the free 
energy density and extracting the Casimir amplitude due 
to the confinement in the time direction. We conclude the 
paper by giving a brief account of our results in the last 
Section. 

Quantum critical behavior. — To gain insights into 
the thermodynamics of model (2) it is more convenient to 
use the field theoretic formulation with the Lagrangian 
density [13-15] 



Jf(^,V) = *t/.^t/' + |^(^tA)^ 



(3) 



The partition function is computed from the path inte- 
gral 



Qixjj&xp exp 
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(4) 



where the fermionic fields are antiperiodic in the imagi- 
nary time direction with period /3 — T^^ i.e. xp{x,l3) — 
~'ip{x,0). Here and below we work in units such as 
h — ks — ^- The discrete chiral symmetry corresponds 
to 

i' = l5ip, ■0 = -'075, (5) 

defining the fermionic condensate {ipip), with thermody- 
namic average taken with the action in eq. (4), as the 
order parameter or mass term of the theory. Under the 
above transformation the Gross-Neveu Lagrangian is in- 
variant. 
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In the large N limit, instead of (3), using a stan- 
dard decoupling technique [15] based on the Hubbard- 
Stratonovich transformation we get the equivalent effec- 
tive Lagrangian: 



(6) 



where A = \{x,t) is an auxiliary scalar field. Let us note 
that in ref. [16] it was shown that the model is renormal- 
isable order by order in expansion. 

The functional integral (6) entering the expression of the 
partition function is Gaussian in fermionic fields. These 
are integrated out to yield the effective action 
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With the aid of eq. (7) we may compute the partition 
function using the saddle point method. Thus we obtain 
an expression for the free energy to the leading order in 
expansion as: 
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where w„ = {2n + 1)'kT are the Matsubara frequencies 
for fermions and ff = {X{x,t)) is the solution of the gap 
equation 
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emanating from equating the derivative, with respect to cr, 
of the free energy j^f{g,T) to 0. For our further analysis 
we will regularise the theory by introducing a cutoff A over 
the spatial wave vector q. 

Eq. (9) has a trivial solution cr = and a nontrivial one 
cr ^ 0. The stability of each solution is to be determined 
by comparing the corresponding free energies. Thus one 
can construct the phase diagram of the model under con- 
sideration. 

To extract the quantum critical behavior of the model 
(6) at zero temperature, we write (9) in a more tractable 
form, employing the Schwinger representation in conjunc- 
tion with the identity 
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that can be derived via the Poisson summation formula 
[17]. Then we have 
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where we have omitted the trivial solution a 
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Setting T = in (11), we find that the quantum critical 
point takes place at the critical coupling gc, to be deter- 
mined from 

1^1 (12) 
8l 2 J (27r)2|q|- ^ ' 

The quantum critical region around gc is defined by the 
condition |g — gcj 1- From the behavior of the lo- 
cal zero temperature fermionic pair correlation function 
{ip{q)ip{—q)) ~ (cr^-l-q^)"^ at the quantum critical point, 
we are able to determine the critical exponent 77 = and 
to identify cr^^ to the correlation length to the leading 
order in the large N expansion. On the other hand the 
mass term, m = (ipip) , is given by 



m 



(13) 



Hence the critical exponent l3 — 1 and consequently v — I 
to the leading order in jj. Other critical exponents can 
be deduced via the hyperscaling relation. 

Now we will turn our attention to the properties of the 
model in a close vicinity of the quantum critical point gc-, 
keeping the temperature finite. For this purpose in lieu of 
eq. (11), we will study the behavior of the mass term at 
finite temperature, say m^, using equation 



-t: ^ = — In 2 cosh — 

£ j?2 2tt 2T 
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obtained via dimensional regularisation. The solution of 
this equation is well known 



, 1 , 

Ut = 2Tarccosh -e 



where we have introduced the scaling variable 
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Let us note that depending on >^, cry can take real or 
complex values. 

Computing the free energy with both the trivial solu- 
tion, cry = 0, and the nontrivial one from eq. (15), we 
find that the phase with the nontrivial solution is stable 
at any T > 0. The free energy corresponding to this phase 
is shown in fig. 1. It has a maximum at T*, which depends 
on the coupling g, where it coincides with that evaluated 
using the triavial solution. In the strong coupling region, 
g > gc, subtracting eq. (11) from its zero temperature 
counterpart we obtain 



m = 2rin2cosh- 



2T 



By setting cr = we get 



m 



21n2' 



(17) 



wherefrom we conclude that the crossover temperature T* 
is proportional to the zero temperature mass term m. We 
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Fig. 1: The free energy as a function of the temperature and 
the distance from the quantum critical point. The thick line 
corresponds to the zero temperature ordered phase, while T* 
shows a crossover line obtained by setting CTt = in eq. (14). 
The crossover line —T* is shown as well. 

would like to notice that T* coincides with the critical 
temperature of ref. [18], where the solution ctt in (15) for 
the variational parameter was assumed to be positive i.e. 
ut > 0. 

Now we will explore the temperature dependence of the 
mass term in the different regions of the {g, T) plane, 
including complex values for Ot as well. In the region 
X — oo, by expanding the r.h..s of eq. (15) we get 

cTT=m-2re-'"/^ + o(e-2m/Ty ^8) 

This shows that in this region we have essentially the zero 
temperature behavior with exponentially small correction 
in temperature. As we approach the quantum critical 
point a crossover occurs at k* = ^ (the T* line on fig 
1). For K < }i* the solution ffy is complex and seems to 
have no direct physical meaning. In particular so long as 
g approaches its critical value gc, the solution Ot creeps 
towards 

^l = -\^'T\ (19) 

This solution shows that the relation between the 
fermionic correlation length and (Jt is not straightforward 
as in the case of zero temperature. To get a meaningful 
result we need to use other means, such as the asymp- 
totic behavior of the local fermionic pair function. Using 
standard techniques [19] we end up with 

« . , I , (20) 

This suggests that the thermal mass term (inverse corre- 
lation length) is defined through 

ct| + TT^T^ (21) 

i.e. the temperature dependent variational parameter is 
shifted by the lowest mode wp- This is in agreement with 



the fact that fermions have no zero mode at finite tem- 
perature [15]. Note that in the limit of zero tempera- 
ture m — limT_>o iTiT coincides with cr = limT~>o ctt as 
expected. A similar situation is encountered in the frame- 
work of the spherical model confined in a film geometry 
under antiperiodic boundary conditions [8,20]. 

The leading behavior of the mass term at finite temper- 
ature in the region >f — > 0^ is 

ruT = ^7rr+ -^7r>fT-HO(>^2). (22) 

By setting g = gc 'we obtain for the critical amplitude of 
the mass term at the quantum critical point the universal 
number -^tt. This shows that the finite temperature mass 
term vanishes linearly in T as we approach the quantum 
critical point from above on the line g = gc vi^ the phase 
diagram. Whence the dynamic critical exponent z for the 
model under consideration is z = 1. 

In the limit x — > — oo, we have Ot ~ i^^T — iTe^'^'^ 
showing that mr vanishes exponentially in the weak cou- 
pling region, g < gc- The behavior of xriT/T as a function 
of >f is shown in fig. 2. It is seen that the universal scaling 
function associated with mx is a monotonically increasing 
function of the scaling variable >f. 
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Fig. 2: Dependence of the scaling functions ^tU<) 2)(>f) 
associated, respectively, to the finite temperature mass term 
(solid line) and the singular part of the free energy (dashed line) 
upon the scaling variable x in vicinity of the critical coupling. 

To sum up, the mass term (21) normalized to the tem- 
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perature is found to behave as: 
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H — > — OO, 

>f = 0, 

>c — )■ OO, 



as a function of the scaling variable >c. Thus we can define 
three different regions in the T) plane depending on the 
behavior of rriT as a function of yc. These are: 

- Renormalized classical: The mass term diverges ex- 
ponentially as we approach the zero temperature phase. 
For T > 0, thermal fluctuations destroy long-range order 
at any finite temperature. 

- Quantum critical: The leading order in T of the mass 
term is linear. 

- Quantum disordered: The mass term is independent 
of the temperature. 

These regions are separated by the crossover lines T* 
and —T* on fig 1. In fig. 3 we show the qualitative 
crossover diagram of model (2). It is worth noticing that 
similar behaviors were found in studies involving many 
other quantum systems (see e.g. ref. [3] and references 
therein) . 
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Fig. 

quantum critical region. The thick line shows the zero temper- 
ature ordered phase, while the dashed lines are the crossover 
lines of fig 1. 

Scaling behavior of the free energy. In the 

quantum critical region, |§ — ^ Ij using the repre- 
sentation 

dt 
T 



ie-'-e'^') 



z > 0, 



(23) 



and dimensional regularisation, the singular part of the 
well known expression for the free energy (8) can be writ- 
ten in the scaling form 



l/.(g,T)=T3 2}(,^), 



(24) 



where we have introduced the universal scaling function 
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in order to make contact with eq. (i). Here cry = (Jt{>c) 
is the solution (15) of the gap equation (9) and 
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are polylogarithms [21]. 

Let us now explore the behavior of the scaling function 
2)(>^). Away from the quantum critical point the leading 
behavior is given by 



i[.2>.+ iC(3)] 



>c — > — OO, 
>c — >■ OO. 



This shows that the scaling function is linear in >c in the 
limit >f — > — OO and decreases as —n^ in the opposite limit. 

In fig. 2 we show the behavior of the universal scaling 
function 2)(>?). It has a maximum estimated to be 



2)(x* 



fsi8,T* 
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C(3) 



-0.143485.. 



a.t m: = m:* = = 0.110318..., corresponding to the 
crossover line T* of fig. 1. 

At the quantum critical point gc we may obtain an an- 
alytic expression using the relations [21] 
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where Clausen's function [22] is defined through 

ch{e) 



(26) 

< 61 < 27r, 
(27) 



sm{n9) 

n=l 



Thus the critical amplitude of the singular part of the free 
energy of the massless Gross-Neveu model is given by 



= -0.274543- •• . 



(28) 



It is negative. In terms of the Casimir effect this implies 
an attractive Casimir "force", despite the fact that the 
boundary conditions imposed over the imaginary time are 
antiperiodic. A similar attractive behavior for the Casimir 
force was found in reference [23,24] induced by the vacuum 
fluctuations of fermionic fields between two parallel plates. 
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Discussion. — We have investigated the quantum 
critical behavior of the massless Gross-Neveu model, as 
a prototype for four-fermionic models. At infinite N it 
has a property that is rarely encountered in the realm of 
critical phenomena, namely exact solvability. This prop- 
erty allows the derivation of thermodynamic functions in 
closed form, without recourse to perturbative methods. 
We have investigated the critical behavior of the model in 
the vicinity of the quantum critical point (zero tempera- 
ture critical coupling) and the changes caused by switching 
on the temperature. In particular we have paid atten- 
tion to the interpretation of the quantum critical scaling 
in terms of the thermodynamic Casimir effect, resulting 
from the confinement in the imaginary time direction. 

Notice that the present investigation rules out any phase 
transition at finite temperature. Such a transition was ob- 
tained in previous studies assuming that the solution (15) 
of the gap equation obeyed by the variational parameter 
of the effective theory can take only positive real values. 
Here such a constraint is released and we consider complex 
values as well. This is justified by the fact that the free 
energy is always real as shown in figs. 1 and 2. In this 
case it was found appropriate to redefine the mass term 
through shifting the variational parameter by the lowest 
fermionic mode. 

The behavior of the free energy and the mass term in 
the vicinity of the quantum critical point was investigated 
in details. We have identified the three different regions 
typical in quantum critical phenomena: Quantum disor- 
dered, quantum critical and renormalized classical (see fig. 
3). These are separated by crossover lines corresponding 
to \g — gc\ ^ T. In particular, we would like to empha- 
size that it was found the mass term is linear in T in the 
region y |g — ~^ 0. The critical amplitude of the mass 
term, which coincides with the inverse of the correlation 
length describing fermionic correlations, was determined 
to be ^TT. 

The temperature dependence of the singular part of the 
free energy in the vicinity of the quantum critical point is 
investigated in details. Its associated scaling function is 
found to exhibit a maximum at The corresponding 
critical amplitude may be interpreted as the Casimir am- 
plitude, due to the confinement in the time direction. This 
is computed exactly [see eq. (28)] and found to be negative 
as it is the case for the quantum nonlinear sigma model 
[25], although this last model obeys the Bose-Einstein 
statistics. This shows that the confinement in the imagi- 
nary time direction is reflected in different way on the sign 
of the Casimir amplitude in comparison with the confine- 
ment in a space direction. 
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